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In order to help detecting superfluidity, we theoretically investigate p-wave pairing super- 
fluids in neutral Fermion atom gases confined by a three dimcnsimentional (3D) harmonic 
potential. The Ginzburg-Landau framework, which is generic for p-wave superfluids, is used 
to describe the order parameter spatial structure, or texture characterized by the ^-vector 
both at rest and under rotation. The Z-vector configuration is strongly contrained by the 
boundary condition due to a trap. It is found that the ground state textures exhibit sponta- 
neous supercurrent at rest both cigar and pancake shape traps. The current direction depends 
on the trapping shape. Under rotation a pair of half-quantum vortex with half- winding num- 
ber enters a system and is stabilized for both trap geometries. We give detailed explanation 
for their 3D structure. The deformations of the condensate shape are seen with increasing 
the rotation speed, which is tightly connected with the underlying vortex formation where 
the condensates are depressed in the vortex core. 
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1. Introduction 

Superfluids with multi-component order parameter (OP) are omni-present, exemplified by 
spinor Bose-Einstein condensates (BEC) 1-4 ** in cold Bosonic atom gases with 23 Na or 87 Rb, 
strongly interacting Fermionic liquid 3 He atoms, 5 ~ 9 ^ and certain heavy Fermion supercon- 
ductors, such as UPt3 10,11 ' 1 or color superconductivity in dense quark- gluon plasmas in high 
energy physics. 12 ' This branch of physics is of great interest from a fundamental physics 
view point as it is expected to exhibit rich topological defect structures or vortices and allow 
investigation of a new phase of matter. Those vortices that can accommodate the Majorana 
zero mode at the core may be useful in quantum computing. 14 ) We will show in this paper 
that half-quantum vortices (HQV) are stabilized, which is a candidate to lead to the Majo- 
rana particles in the core. 14 ) Therefore the present p-wave superfluids in general provide a 
tesing ground to explore rich physics associated with topological defects. Previously we have 
demonstrated that superfluid 3 He-A phase is another fertile ground to find the Majorana par- 
ticle. 15, 16 ) This is the second example to show HQV as a stable vortex in our study 16 -* where 
superfluid 3 He-A phase confined between parallel plates is demonstrated to exhibit HQV. 



*E-mail address: tsutsumi@mp.okayama-u.ac.jp 



1/18 



J. Phys. Soc. Jpn. Full Paper 

Recently, p-w&ve resonance superfluidity attracts much attention in Fermionic alkaline 
atom gases, such as 6 Li 17 ~ 19 ) and 40 K, 20 ~ 22 ) both experimentally and theoretically. 23-27 ) Ex- 
periments into achieving p-w&ve resonance superfluidity are steadily progressing 19 ) and hence 
it is timely and necessary to consider the generic properties of p-w&ve superfluidity both at 
rest and under rotation to detect its superfluidity, which is signified by non-trivial topological 
structures or vortices. In this respect, the interest is in the spatial structure, i.e. texture, of 
the nine component OPs describing the spin triplet p-w&ve superfluidity. A p-w&ve Feshbach 
resonance occurs at the different magnetic field in each hyperfine spin state of Cooper pair. 
Since the spin state of superfluidity is fixed by the external magnetic field, the spin degrees 
of freedom are frozen, and hence only the orbital degrees of freedom are active. The order 
parameter space consisting of orbital three components is analogous to that of superfluid 3 He, 
in particular, the A phase to which spin and orbital state can be separated, 5 ~ 8 ) where the OP 
is described by a tensor 

= d^Ai = x,y,z). (1) 

The and Ai describe the spin and orbital state of a Cooper pair, respectively. Here the OP 
is characterized only by Ai, described as 

A(p) = A xPx + A y py + A z p z . (2) 

In a sense, a p-w&ve superfluid is analogous to the "spinless" superfluid 3 He-A phase. 15 -* 

The dipole-dipole interaction between two alkaline atoms acts to split the relative orbital 
state for two particles, depending on the projections of the orbital angular momentum, either 
mi = il or mj = 0. This results in breaking of the degeneracy between 

p± = ^-^={p x ± ipy) and p = p z . (3) 

This splitting was estimated to be large for 40 K by Cheng and Yip, 23 ) evidenced by the clear 
difference in the Feshbach resonance mag netic fields (splitting field = 0.47±0.08G). 21 ) For 6 Li, 
the splitting may be small, as an experiment conducted in a magnetic field of H = 158.5(7)G 
shows no clear resonance splitting. 19 ) 

A critical difference between superfluid 3 He and a p-wave resonance superfluid of atom 
gases lies in the boundary conditions. In the superfluid 3 He-A phase the /-vector, which is 
the orbital angular momentum of Cooper pairs, is always perpendicular to a hard wall so 
that the perpendicular particle motion is suppressed. In other words, the point nodes in the 
/-vector direction touch the hard wall so as to minimize the condensation energy loss at 
the boundary. 6 -* On the other hand, atom gases are confined by a three dimensional (3D) 
harmonic trap potential, where the condensation energy density gradually decreases towards 
the outer region, the Z-vector tends to align parallel to the circumference. This orientation 
is advantageous because the condensation energy is maximally gained by allowing the point 
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nodes to move out from the system. The trap potential is easily controlled, resulting in various 
shapes, such as cigar or pancake shapes. As we describe below, the trapping potential can be 
an important tool to control the 3D texture. Indeed, the 3D trapping structure constrains the 
possible textures. Our purpose is to investigate a possible 3D textures in a 3D harmonic trap 
potential and thereby to help identifying p-wave superfluidity. 

The organization of this paper is as follows: We employ the Ginzburg-Landau (GL) frame- 
work which relies only on global symmetry principle. The GL free energy functional form is 
introduced and the relevant physical quantities, such as supercurrent and ^-vector are given 
in § 2. We show the phase diagram of the stable states for uniform and infinite system, and 
explain how to numerically examine realistic confined systems in § 3. Section 4 presents the 
stable texture in the cigar shape trap, especially, we mention the spontaneous supercurrent at 
rest and the half-quantum vortex (HQV) under rotation. In § 5, we show the different textures 
stabilized in the pancake shape. The final § 6 is devoted to a summary and discussion. A short 
version of the present paper is to be published. 2S ) 

2. Formulation 

Here we employ the GL framework. 6 ) This framework is general and flexible enough to 
allow us to examine a generic topological structure, and applicable to cold Fermionic atom 
gases with a harmonic trap potential under ksT c S> huj, where T c and to are the transition 
temperature and the trap frequency, respectively. 29 ^ In terms of the tensor A^i forming OP of 
p-wave pairing the most general GL functional density /bulk f° r the bulk condensation energy 
up to fourth order is described as 

/bulk = —otiA^A^i + fixA^A^AyjAyj + f32A fMi A /1 iA U jA U j+f3sA fli A ui A^jA l/ j 

+(34 : A jli A V iA u: jA^j+j3^A fli A U iA U jA jl j, (4) 

which is invariant under spin and real space rotations in addition to the gauge invariance 
U(l)xSO^(3)xSO ( ' L )(3). The fourth order terms are characterized by five independent in- 
variants, (5\ ~ /?5 in general. 6 ) Since the spin degrees of freedom are frozen due to applied 
magnetic field for magnetic Feshbach resonance, only the orbital degrees of freedom Ai in 
A^i = d^Ai are active. Namely it reduces to 

/bulk = -ao(l - U)A*Ai + (3 2A A*A t A*A 3 + ^A*A*AjA 3 , (5) 

where /?24 = P2+P4 and ti = T/T c i (T c j is the transition temperature for the i-component). As 
mentioned, the dipole-dipole interaction causes splitting of the transition temperatures into 
two groups T cx = T cy and T cz . We introduce a = T cx /T cz , which indicates the degree of the 
broken symmetry of the system and characterizes atomic species used. 21 ) The pairing state 
having the orbital projection mi = is favorable over mi = ±1, namely < a < 1 due to 
the dipole-dipole interaction. 21 ) The three components become degenerate for a — > 1. When 
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a — > 0, the polar state with the OP A(p) = A z p z tends to be stable. 

The gradient energy consisting of the three independent terms 6 -* is given by 

/ gr ad = *i WMJXftA;) + KtWAZWjAi) + K,(dtA\){d 3 A 3 ). (6) 

The centrifugal potential energy due to rotation with fi, which is derived in Appendix A, is 
written as 

2 

/cent = -^n 2 p 2 (K 1 A*A i + K 2 \A e \ 2 + K 3 \A e \ 2 ), (7) 

where d{ = V, — i(m/h)(ft x r)j. For fi || z, p 2 = x 2 + y 2 and ^ = — AjSinfl + A^cosfl in 
the cylindrical coordinates. 

The GL parameters ao, /?24 = 02 + Pi, ft and K\ = K 2 = K$ = K are estimated by 
taking the weak coupling approximation, assuming the Fermi sphere: 6 ^ 

N M R R R 7 C(3)iV(0) _ a m 

ao = ^-' ft = ft = ft = i 2 o (7rfcB r c) 2 = P (8) 

and 

7 C(3)iV(0)(fa F ) 2 
240(vrA; B T c ) 2 

where iV(0) is the density of states at the Fermi level and vf is the Fermi velocity. The weak 
coupling approximation should be a good guide for understanding the generic properties of 
the p-wave superfluids of atom gases because it has been applied successfully, even to liquid 
3 He with strong interacting Fermions and only small additional strong corrections. 6 ) 
It is convenient to discuss (5), (6) and (7) in the following dimensionless units, 

^*'ih r 'i^ n (10) 
with the zero-temperature GL coherence length 



k = ;t 2 ( q ) 



The units of OP and the angular velocity are 



« m ~ iw <13) 



/bulk = - (1 - U)A*Ai + A*AiA*A 3 - + -A*A*A 3 A 3 , (14) 



In the dimensionless unit, (5), (6) and (7) are written as 

1 

2" 

/ grad =(d*A*)(d t A 3 ) + (d*A*)(d 3 A t ) + (d;A*)(d 3 A 3 ), (15) 
f ccnt =-n 2 p\A*A i + 2\Ae\ 2 ), (16) 
respectively, where d% = Vj — x r)j. From now on we use the dimensionless expressions. 
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The harmonic trap potential term 29 ) is 

/harmonic = LU 2 ± (p 2 + \ 2 Z 2 )A* A h (17) 

where the dimensionless radial confining potential is uj± and the anisotropy of the harmonic 
trap is expressed as A = ll> z /ujj_. The harmonic trap potential term acts to lower the transition 
temperatures. It is interesting to note that the centrifugal potential leads to the non-trivial 
form, because the OP label implies the orbital angular momentum, a feature absent in a 
spinor BEC. 30 ) The extra factor of 2 1 | 2 in the above form of (16) becomes important when 
evaluating the critical angular velocity Q CI , above which the superfluid flies apart. That is, 
r2 cr = u±/y/3 is greatly reduced from the usual case (Q, cr = u>±). 

The relevant physical quantities are described in terms of the OP Ai(r) as follows. The 
total free energy: 

F = J d ? 'r (/bulk + /grad + /cent + /harmonic) • (18) 

ji(r) = 21m [A)ViAj + A*VjAi + A*V 3 A 3 ] . (19) 



The current density: 



The ^-vector: 



with 



= -^i*r^yj2 ( 2 °) 



| A(r) \ 2 = A* A. (21) 

3. Preriminary Considerations 



Before considering a realistic confined system, we first investigate an infinite system at rest. 
Then the state minimizing the condensation energy (14) is realized. As shown in Appendix B, 
the phase diagram in Fig. 1 consists of the three phases, A, B, and normal (N) phases. The B 
phase is described by A(p) = A z p z , i.e. the polar state. The A phase is described by a chiral 
OP expressed by 

A(p) = A z {p z +i 1 p A _), (22) 

where (0 < 7 < 1) with j5j_ = p x cos<p + p y sm.(j) (4> is arbitrary). This phase breaks the time 
reversal symmetry. The value of 7 depends on temperature and anisotropy a, 

= 2 - (3 - a)t x 
7 2 + (l-3a)t x " 1 j 

In the absence of the dipole-dipole interaction (a = 1), the value of 7 is unity. As the dipole- 
dipole interaction increases (a decreases), 7 decreases. The decrease of 7 becomes larger at 
high temperature. When the value of 7 vanishes, the second order transition from the A phase 
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Fig. 1. (Color online) Phase diagram of the p-wavc pairing state in an infinite system at rest, showing 
temperature (t x — t y ) versus anisotropy a = T cx /T cz . N: normal state, A: chiral state A z (p z +i"/pj_) 
and B: polar state A z p z . 

to the B phase takes place. The transition temperature t c is given by 



The phase diagram corresponds to the case being not so large Feshbach resonance splitting in 
BCS regime. 24 ^ However, it is warned that the determined phase boundary is of qualitative 
at low temperatures which is beyond the GL framework. In the following we examine the A 
phase in confined geometries. 

In order to obtain stable texture of the condensates in a realistic harmonic potential, 
we have identified stationary solutions by numerically solving the variational equations: 
5f(r)/5Ai(r) = in three dimensions where /(r) is the GL energy density functional, the 
integrand of (18). We start with various initial configurations, including singular vortex state 
and non-singular vortex state, and determine the most stable texture by comparing the total 
GL energy (18). 

4. Cigar Shape Trap 

We first consider the stable texture for a cigar shape trap with the trap anisotropy A = 0.2. 
We take 80 x 80 x 120 meshes with the cloud sizes for the Thomas-Fermi approximation 
Rx = Ry = 10 and R z = 50. We fix the temperature at t x = t y = 0.4 and the anisotropy 
parameter a = 0.9. 

4-1 Stable texture at rest 

The stable /-vector texture at rest is shown in Fig. 2. Figure 2(a) displays the amplitude 
distribution of the /-vectors. It can be seen that the amplitude |1| is maximum in the central 
region, and towards the outer regions |1| decreases gradually. At the top and bottom ends, 
the polar state is realized where the /-vector vanishes. Three cross sections are shown in Figs. 
2(b)-(d). In Fig. 2(c), which corresponds to the middle cross section, the /-vectors lie in the 
x-y plane, showing a streamline type pattern in which the /-vectors follow the circumference, 
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Fig. 2. (Color online) Stable texture at rest for cigar trap A = 0.2. (a) Distribution of |1| in the z-x 
plane, (b)-(d) Three cross sections indicated in (a), showing l z (color) and l x and l y components 
(arrows), where the dot and cross marks indicate the imaginative source and sink for the Z-vector 
stream lines. 

like a fluid streaming along a circular boundary. Outside the condensates, an unseen sink 
and source of the Z-vectors exist, giving two imaginary focal points situated outside. Namely 
the left dot and right cross marks in Figs. 2(b)-(d) correspond to the source and sink where 
the /-vectors appear and disappear. This streamline like texture contrasts with the so-called 
Pan- Am texture in superfluid 3 He-A phase 6 ) where the /-vectors tend to point perpendicular 
to the wall due to the boundary condition. In the upper (Fig. 2(b)) and lower (Fig. 2(d)) cross 
sections the streamline texture is maintained, but an l z component appears in addition. 

The associated supercurrent structure is depicted in Fig. 3. The j z component shows 
a circulation supercurrent along the z axis (see Fig. 3(a)). Since f2 = 0, this circulation 
supercurrent is spontaneously generated. This non-trivial condensates flow can be explained 
in supercurrent characteristic of chiral p-wave superfluid: 

j = p s v s + C(V x 1), (25) 

where p s and C are diagonal tensor coefficients. The supercurrent has a structure that is 
similar to that for the total charge-current density of a system in classical electrodynamics. 
The curent in classical electrodynamics is composed of two parts, the actual charge transport 
and an effective current in proportion to V x M due to the magnetization M generated by the 
internal motion of the electrons. In chiral p-wave superfluid the Cooper pair may be thought 
of as representing the atom and 1 the magnetic orbital moment M. Therefore the first term in 
the supercurrent describes the usual flow of Cooper pairs, while the second one is an "orbital" 
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Fig. 3. (Color online) (a) Spontaneous circulating current flows at rest (f2 = 0) in the z-y plane along 
the z direction. Three cross sections indicated in the left figure, showing j z (color) and j x and j y 
components (arrows). (b)The current pattern under rotation (f2 = 0.3). It is seen that the in-plane 
components j x and j y increase, producing the in-plane circular current due to rotation. 



supercurrent driven by the spatial variation of the /-vector. 

In this system, the spontaneous supercurrent is not a usual supercurrent p s v s but an 
"orbital" supercurrent C(V x 1) with the trap potential. In the middle cross section in Fig. 
3(a), the /-vector in-plane bending (V x \) z produces a perpendicular current j z . However, at 
the upper and lower planes in Fig. 3(a), the supercurrent acquires the j x and j y components 
because of the non-vanishing l z component. Therefore, the perpendicular current at the center 
bends such that the condensates are conserved. It is clear from Fig. 3(a) that the supercurrent 
circulates perpetually along the z direction parallel to the long axis of the trap. This result is 
non-trivial and a remarkable manifestation in the topological nature of the texture. 

In Fig. 3(b) we also display the current patterns under rotation (f2 = 0.3u;_l) for com- 
parison at rest. It is seen that the in-plane components j x and j y increase, producing the 
in-plane circular current due to rotation, in particular in the middle cross section. Thus under 
rotations, the current consists of the spontaneous one along the z direction and the induced 
circular current in the plane. 

4-2 Half- quantum vortex under rotation 

In Fig. 4 where we depict the middle cross sections of the stable solutions for various 
rotation speeds (also see Fig. 6). It is seen that under rotation, the 3D texture deforms 
continuously and smoothly. As the rotational speed increases, the /-vectors in the x-y plane 
pointing in the x direction acquire a negative z component, as seen by the color change from 
green to blue. This deformation yields in plane circular "orbital" supercurrent as already 
shown in Fig. 3(b). 
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Fig. 4. (Color online) (a) Texture change with rotational speed f2 relative to the in-plane trap fre- 
quency cjj_. Middle cross sections are displayed showing l z (color) and l x and l y components 
(arrows). As 51 increases, I- vectors acquire a negative z component (seen as a color change). At 
Q = 0Alu± two plus HQVs enter from the y direction seen as yellow objects. At 51 = 0.5cl>j^ four 
plus HQVs and two minus HQVs are present and the condensates expand and deform, (b) The 
enlarged figure of HQV where the local coordinates (r, 9) are shown. 



Above a certain rotational speed (f2 = 0Au>±), two HQVs enter from the y direction, where 
the ^-vector at the core pointing in the positive z direction, as seen by the yellow objects. The 
OP far away from the HQV core is described in terms of the local coordinated (r, 9) centered 
at the core (see Fig. 4(b)) as 

A(r » 1,0, p) = exp[i(6/2 + Tr/2)]\A xy \[ S m(6/2)p x - cos(0/2)p y ] + \A z \p z , (26) 

where \A xy \ and \A Z \ are amplitude of the polar state characterized by the basis functions 
Px-Py [x-y polar state) and by p z (z polar state), respectively. The HQV is formed by only 
the x-y polar state as is seen from (26). The above state (26) can be written as 

A(r» 1,9, P ) = -^=\A xy \{p + + e i6 p„) + \A z \p . (27) 

The vortex core of this HQV has the orbital angular momentum with positive z component 
because p + is non-vanishing there, so we define it as a plus HQV. Similarly, one can construct 
a minus HQV with the identical vorticity but the local orbital angular momentum of the 
negative z component in its core (p—)- 
The realized OP written as 

A(r, p) = A + (r)p+ + A_(r)p_ + A (r) Po , (28) 
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A ± = T^{A x TiA y ),A = A z (29) 

is shown in Fig. 5(a) where the amplitudes (upper line) and phases (lower line) for each 
component A + ,A^,Aq are displayed for 17 = 0Alo± corresponding to Fig. 4. It is seen that 
there are two vortices with the winding number 1 in A^ component on the y axis. The 
depletion of the OP amplitude by the vortices is compensated by the growth of A + amplitude. 
The HQV is embedded in the surrounding OP field. As we walk around one of the HQV core, 
the p-wave pairing state changes in the following manner: At 6 = it, the OP is described 
as A(p) = — \A xy \p x + \A z \p z , namely this is the polar state. At 6 = 0, the OP is now 
A(p) = —i\A xy \py + \A z \p z . This is a superposition of the polar state and the chiral state. The 
spatial position where the arrow of ^-vector near the HQV core vanishes is the polar state, 
and on the opposite side of HQV core, the chiral state is superposed by the polar state. Since 
this HQV breaks the reflection symmetry, the isolated HQV is energetically disfavored and 
rather a pair of the HQV is advantageous. This is one of the reasons why we found a pair of 
the HQV at Q = 0.4oj ± . 

Upon further increase in the rotational speed (Q = 0.5wj_), the HQVs enter further from 
the x direction (see Fig. 4). They are different from above mentioned plus HQVs. As is seen 
from Fig. 5(b), a pair of plus-minus HQVs appears from the x direction. There are the vortices 
with the winding number 1 in each A + and A- component on the different position of the 
x axis. The depletion of the OP amplitude by the vortices is compensated by the growth 
of amplitude with the opposite orbital angular momentum component. The minus HQV is 
situated more inside in the trap potential than plus HQV, because of the repulsive interaction 
between the plus HQV on the x axis and the y axis. It is an analogous situation for the HQVs 
in the F = 1 spinor BEC. 31 ) We also notice from Fig. 4 (f2 = 0.5wj_) that the condensate 
profile itself expands, deforms and deviates clearly from a circular form due to the non-trivial 
centrifugal energy. 

Figure 6 shows a different view of Fig. 4, displaying the z-x cross section. At rest, the 
^-vectors point almost to the x direction. As Q increases, the downward l z component appears, 
which causes a counterclockwise circular "orbital" supercurrent. At the top and bottom ends 
of the system the HQVs appear as indicated by asterisks. These blue lines show the polar 
state neighboring the HQV as mentioned above. In the Q, = 0.5u>± case, the side view of two 
pairs of plus-minus HQVs can be seen clearly. The light blue lines indicated by the red arrows 
correspond to the polar states between plus and minus HQV. These lines bend outward away 
from the center towards ±z direction. Around the positions indicated by the red arrows the 
condensate profile is greatly deformed. Because there the polar state dominates over the chiral 
state, which leads to the depletion of the condensate, resulting in this deformation. 
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Fig. 5. (Color online) The amplitude (upper) and phase (lower) of the OP components A+, A- and 
Ao in the x-y plane at z = corresponding to Fig. 4. (a) At f2 = 0.4cjj_ the vortices with winding 
number 1 enter A- component. The depletion of the OP amplitude by the vortices is compensated 
by the growth of A + amplitude, (b) At £1 = Q.5uj± the vortices with winding number 1 appear. 
The vortices on the y and x axis are plus and plus-minus HQVs, respectively. 
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Fig. 6. (Color online) The z-x cross sections corresponding to Fig. 4. The Z-vectors lie almost in the 
x-y plane at rest. Under rotation two HQVs appear, shown as asterisks. At ft = 0.5wj_ two pairs 
of the HQVs can be seen, indicated by red arrows. 
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5. Pancake Shape Trap 

In order to understand the geometric effect of the textures stabilized in cigar shape trap, 
we examine the pancake trap case. We consider the anisotropy A = 3.0 as an example, where 
we take 100 x 100 x 80 meshes with the cloud size R x = R y = 30 and R z = 10 in the 
Thomas-Fermi approximation . We fix the temperature at t x = t y = 0.5 and a = 0.95. 

5.1 Axis symmetric texture at rest 

Figure 7 displays the resulting /-vector texture (Fig. 7(a)) and associated in-plane super- 
current (Fig. 7(b)) at rest. It is seen that most of the /-vectors point to the negative z direction, 
except for those near the upper and lower surface regions, which acquire the r-component of 
the radial direction. The trap potential forces the /-vector to be parallel to the surface of 
condensates. This effect is especially strong when the curvature of the surface is large. Thus, 
in this case shown in Fig. 7(a), the left and right ends of the system force the /-vectors to 
point to the z direction, giving an overall /-vector configuration to the z direction, even for 
the vectors near the center. The amplitude of the /-vectors decrease, namely the polar state 
mixes with the chiral state towards the outside. 

This stable texture is axis symmetric around the z axis, which is different in the cigar 
shape trap. This is because the effect of the trap potential is greater than one of the dipole- 
dipole interaction in the present situation (a = 0.95), so that the chiral state consisting of the 
p x and p y components dominates over the polar state, resulting in the axis symmetric texture 
with respect to the z axis. If the influence of the dipole-dipole interaction is much greater, or 
a becomes small, the stable texture is similar to those in the cigar shape trap. 

The spatial variation of the amplitude of /-vectors towards the outer region generates a 
circular "orbital" supercurrent in the x-y plane as shown in Fig. 7(b). Since the spontaneous 
supercurrent flows around the external rotational axis, the axial symmetric texture is stable 
against low rotation. 

5.2 Vortices under rotation 

As shown in Fig. 8 where the x-y cross section at z = is displayed for various rotations. 
Similarly to the previous texture changes in the cigar shape, we see a pair of the HQV enter 
from the y axis for fi = 0.3wj_. Those vortices are exactly the same HQV as in the cigar case 
(see Fig. 4). The following sequence upon increasing O is very similar to the previous cigar 
case, slightly differing its rotation speed where in the pancake case the sequence is shifted 
to lower speed. We illustrate the x-z cross section at y = in Fig. 9. The almost /-vectors 
pointing to the negative z direction at rest now deforms as Q increases. From the right and left 
ends the polar state seen as blue color invades into the system. At CI = 0.5wj_ the condensate 
profile itself is deformed greatly because of the dominance of the polar state associated with 
increasing HQ Vs. 
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Fig. 7. (Color online) Stable texture in the pancake trap with A = 3.0 at rest, (a) /-vector pattern in 
the z-r plane. The pattern is axial symmetric around z. (b) Spontaneous supercurrent in the x-y 
plane at 2 = 0. 



6. Summary and Discussions 

By minimizing the generic GL energy functional within weak coupling approximation, we 
find stable textures for p-wave superfluids in neutral atomic gases confined in 3D harmonic trap 
potentials which are to be realized in near future. In order to help detecting its superfluidity, 
we examine two typical trap geometries, cigar and pancake shapes. At rest, the obtained stable 
textures in both shapes exhibit the spontaneous supercurrent flow. Under rotation, a pair of 
the HQVs enters in the condensates. The isolated HQV is never stabilized in our calculations 
because of topological constrains. 

It is interesting to note that the direction of the spontaneously generated supercurrent at 
rest is always perpendicular to the direction of the majority /-vectors, that is, the condensates 
in the cigar and pancake case the /-vectors lie on the x-y plane and point to the z direction, 
respectively, so the supercurrent flows toward the z direction and in the x-y plane. This implies 
that the trap shape is critical in understanding and controlling the physics of the textures on 
the p-wave superfluid. 

The p-wave Feshbach resonance occurs by using either the same species 23 ' 24 ) or two species 
with different hyperfine states 25-27 -* whose difference in our context amounts to giving different 
values of the parameter a because the dipole-dipole interaction for a Cooper pair, yielding 
the splitting between mi = ±1 and mi = 0, works differently for two cases. 

As for Major ana particle generated in the vortex core in half-quantum vortex, let us 
examine whether or not the present HQV accommodates the Majorana particle in its core. 
The Mojorana condition is obviously fulfilled when the half-quantum vortex is involved the 
d-vector, namely the spin degrees of freedom. 16 ' However, in the present spinless HQV it turns 
out by examining the Bogoliubov-de Gennes equation 32 - 1 that the particles bound in the HQV 
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8. (Color online) The x-y cross sectional view of the textures for the pancake shape trap with 
A = 3.0 for various rotations. The two objects on the y axis for f2 = Q.3u>± are the HQVs. Compare 
those with Fig. 4 in cigar shape case. 




9. (Color online) The x-z cross sectional view of the same textures in Fig. 8. The polar state 
indicated by blue color gradually enter from the right and left ends, making the system deform. 
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core always have a finite energy, not at zero-energy, which is indispensable for the Major ana 
condition. Therefore, those are never the Majorana particles. In order to find the Majorana 
particles, it is necessary to stabilize the singular vortex with odd winding number in chiral 
superfluids. 32 ) In the present case, it might be realized for extremely flat pancake limit where 
the p z component becomes irrelevant. The detailed energetics between the present HQV and 
the singular vortex belongs to future study. 
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Appendix A 

We start with the Hamiltonian in a co-rotating frame: 

H = H - n ■ L, (A-l) 
#o = £ [(l>?/2m) + U( ri )] + V (A2) 

i 

where Hq is the Hamiltonian in a non-rotating system, consisting of the kinetic energy term, 
the harmonic trap potential term U and the interaction energy term V. The angular velocity 
due to the external rotation is ft and the angular momentum is L = J2i 1/2 (fj x pi — pi x r,). 
We can write H in the form 

H = E (P* " ™n,i) 2 + U{n)] +V-Y: \m<,i, (A-3) 



2 



where v nj j = ft x r« is the "normal fluid" velocity at the location of the particle i. The last 
term is the centrifugal energy. 

The gradient energy given by the first term of the Hamiltonian (A-3) 

fgrad = 16{kIbTc) 2 K p * ~ mVn ^ A ^ ^ P ~ mVn ^ fc A d ( VFi Pi VFk P l h ' ( A " 4 ) 
where (■■■)$ denotes the Fermi surface average. Taking the contraction of the subscripts in 
the Fermi surface average, the mean value is finite. Replacing p with —ihV, we obtain the 
well-known form of the gradient energy 6 ) 

/grad = Ki(d*Aj)(diAj) + K 2 {d*A*){d 3 A t ) + K^A*)^). (A5) 

In a similar manner, the centrifugal energy given by the last term of the Hamiltonian (A-3) 
can be recast into 

/cent = - iQ^foT) 2 t( mv ")^j] [(mv n ) k Ai] (v F iPjV Fk pi) p (A-6) 

2 

ill 

= --^[K 1 (nx r)iA*(n x r)iAj + K 2 (n x r)iA*(Cl x v) j A i 



15/18 



J. Phys. Soc. Jpn. Full Paper 

+ K 3 (ft x r)iA*(ft x r)jAj}. (A-7) 
Taking fi || z, and 12 x r = we finally obtain the expression of the centrifugal potential 



as 

2 

/cent = -^Vt 2 p\K x A\A t + K 2 \A e \ 2 + K 3 \A e \ 2 ). (A-8) 

Appendix B 

The bulk free energy (5) is written as 

/bulk = - (1 - t x )(\A x \ 2 + |A/) - (1 - at x )|A*| 2 

+ (IA,! 2 + \A y \ 2 + |A 2 | 2 ) 2 + l -(Af + Af + Af){A 2 x + A 2 + A 2 ). (B-l) 
By taking A x real, A y = \A y \e l9y and A 2 = \A z \e t9z without loss of generality, we write it as 
/bulk = - (1 - t x ){\A x \ 2 + \A y \ 2 ) - (1 - at x )\A z \ 2 

+ ^(|A X | 4 + \A y \ 4 + |A 2 | 4 ) + 2(\A x \ 2 \A y \ 2 + |A/|AJ 2 + lAJ 2 ^^ 2 ) 

+ \A x \ 2 \A y \ 2 cos29 y + \A y \ 2 \A z \ 2 cos(29 y - 29 z ) + \A Z \ 2 \A X \ 2 cos26 z . (B-2) 

There are six cases, which are possibly the minimum solutions: (i) 9 y = 9 Z = 0, (ii) 9 y = 0, 
9 Z = tt/2, (iii) 9 y = vr/2, # 2 = 0, (iv) 9 y = vr/2, 9 Z = vr/2. (v) 9 y = vr/3, 6> 2 = 2vr/3 and (vi) 
|A X | 2 = |A^| 2 = 0. 

The relevant minimum solutions are found for (ii) and (vi). In the former case (ii) the 
solution is given by 

/bulk = -^[3(1 " tx? ~ 2(1 - t x )(l - at x ) + 3(1 - at x ) 2 ] 

1 2 i \ a 1 2 1 ■ 



|A^ + \A y \ z =-[3(1 - «*) - (1 - at x )] (B-3) 
|A 2 | 2 =^[-(1-^+3(1 -at,)] 



This solution is valid for 2 + (—3 + a)t x > 0, which determines the boundary between this 
phase called the A phase and the single component phase called B phase below. The other 
solution for (vi) is expressed as 

/bulk = _ T^(! - at x) 2 

l \ (B-4) 
|A 2 | 2 = -(l-at x ) 

This B phase is described by a single component, thus it corresponds to the so-called polar 
phase. 
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